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The Lagrangian point of view is adopted to study turbulent premixed combustion. The evolution 
of the volume fraction of combustion products is established by the Reynolds transport theorem. It 
emerges that the burned-mass fraction is led by the turbulent particle motion, by the flame front 
velocity, and by the mean curvature of the flame front. A physical requirement connecting particle 
turbulent dispersion and flame front velocity is obtained from equating the expansion rates of the 
flame front progression and of the unburned particles spread. The resulting description compares 
favorably with experimental data. In the case of a zero-curvature flame, with a non-Markovian 
parabolic model for turbulent dispersion, the formulation yields the Zimont equation extended to 
all elapsed times and fully determined by turbulence characteristics. The exact solution of the 
extended Zimont equation is calculated and analyzed to bring out different regimes. 
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Turbulent premixed combustion is a challenging scien- 
tific field involving noncquilibrium phenomena and play- 
ing the main role in important industrial issues such as 
energy production and engine design. 

A Lagrangian point of view is here adopted, lead- 
ing to a description of turbulent premixed combustion 
which takes into account, for all elapsed times, the tur- 
bulent dispersion, the volume consumption rate of re- 
actants, and the flame mean curvature. The proposed 
approach generalizes and unifies classical literature ap- 
proaches that are based on the so-called level-set method 
PQ or are based on the Zimont balance equation, origi- 
nally hinted at by Prudnikov [2] , also known as Turbulent 
Flame Closure model [3J . Moreover, the proposed formu- 
lation has the striking property to be compatible with 
every type of geometry and flow in an easier and more 
versatile way than previous approaches, and it emerges to 
be easily modifiable to include more detailed and correct 
physics. It is worth recalling that the Zimont equation 
was introduced on the basis of experimental observations 
and that a great effort has been undertaken to give a 
deeper theoretical foundation to it [3-6J. The present 
Lagrangian formulation constitutes a reliable theoretical 
support for the Zimont combustion model. 

The process of turbulent premixed combustion is 
mainly characterized by flame propagation towards the 
unburned region and turbulent dispersion of the resultant 
product particles. The combustion process is described 
by a single dimensionless scalar observable, denoted as 
average progress variable, < c(x,t) < 1, and represent- 
ing the burned-mass fraction, i.e., the fraction of burned 
particles which are located in x at time t. The value 
c(x, t) = 1 describes the presence of only products and 
the value c(x, t) = describes the presence of only re- 
actants. To avoid unnecessary mathematical difficulties, 
we consider a constant-density mixture and a zero-mean 
turbulent velocity field. Molecular diffusion is also ne- 
glected. 

In this Letter, the fresh mixture is intended to be a 
population of particles in turbulent motion that, in a 



statistical sense, change from reactant to product when 
their average positions are hit by the flame. Let Vl{t) 
be the portion of space surrounded by the flame surface; 
then those particles with average position (x) G il(t) 
are marked as burned particles. The occurrence in x at 
time t of a particle transit is described by a probability 
density function (PDF). Let p(x; i|xo) be the PDF asso- 
ciated with a particle displacement where Xo is the initial 
condition of a Lagrangian trajectory and, without loss of 
generality, let t = be the ignition instant. With the 
assumption that particle trajectories are not affected by 
the chemical transformation, the average progress vari- 
able c(x, t) turns out to be defined as the superposition 
of PDFs of burned particles, i.e., those p(x;t\xo) with 
(x) e £!(i). For a zero-average velocity field, the particle 
average position is (x) = xo and then 

c(x,i) = / p(x;i|x )ebc . (1) 
Jn(t) 

The evolution law for the progress variable c(x, t) is ob- 
tained applying Reynolds transport theorem to ([I]) which 
gives 

3c f Op f 

^7 dx o + / V *o ' [ u ( x o, 0p( x ; *|xo)] dx , 
at Jn(t) at Jn(t) 

(2) 

where V Xo is the gradient with respect to Xo and u(x, t) 
is the expansion velocity field of ft (t). 

Let the turbulent dispersion be represented by the gen- 
eral evolution equation 

dp 

Oj:=£x[p]> p(x;0|xo) = «y(x-Xo), (3) 

where the spatial operator £ x [-] includes the particle dis- 
placement statistics such as the variance u 2 (t) = (||x — 
xo|| 2 >/3. 

Equation (T2J) is also governed by the volumetric ex- 
pansion of Q(i). This expansion is connected with the 
consumption rate that in a general form is set to be 

Vc 

u(x,i) =U{K,t)h, h = -jT—rT, (4) 
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where k(x, t) = V • h/2 denotes the local mean curva- 
ture. Since molecular processes are neglected, the initial 
burning speed must be zero, i.e. U(k, 0) = 0. From Q 
the location of the flame surface follows to be 



(5) 



£ f (t) = £ + f u(C f ,T)dT. 
Jo 

Finally, inserting ^ and Q in ([2| gives 
£ x [c] + / u V XQ p dx 

JQ(t) 

V Xq k • h + 2U(k, £)k(x , t) dxo ■ (6) 



dc 
df 



n(t) 
dU 



The evolution of the progress variable is then led by three 
factors: turbulent motion, displacement speed of the con- 
tours of c(x, t) and their mean curvature. It is worth re- 
marking that ^ cannot be reduced to the most widely 
used front propagation equations [7j, and, since ^ fol- 
lows from the exact Lagrangian definition none of 
them is physically correct to model turbulent premixed 
combustion. 

When particle motion is neglected, products and re- 
actants turns out to be frozen, dp/dt — 0, so that 
p -> S(x - x ) and J n(t) V Xo • [u(x , t)p}dx = u ■ Vc. 
Here the identity Vp = — V Xo p has been used. In this 
limit case, using Q, Eq. ^ reduces to 



Q { l, ifxen(t), 

-=U{K,t)\\Vc\\, c(x,i) = { (7) 
at { 0, otherwise, 

the celebrated Hamilton-Jacobi equation stated by 
Sethian [T] to track the flame front surrounding the 
burned volume O(t), which is related to the G-equation 
[8] and to the Kardar-Parisi-Zhang equation [9]. Equa- 
tion |7]) can be now interpreted as a consequence of 
Reynolds transport theorem. 

When the normal n to the contours of the progress 
variable is assumed constant, then the mean curvature 
k is zero. Assuming an homogeneous, isotropic and 
stationary turbulence, for the Lagrangian PDF it holds 
Vp = — V Xo p and setting U(0, t) = I4{t) formula ^ turns 
out to be 



dc 
Of 



= £ x [c]+U(t)\\Vt 



(8) 



It must be observed that in ([8]) the turbulent dispersion 
and the flame expansion enter in the progress variable 
evolution with the particle displacement variance u 2 (t) 
and I4(t), that are so far independent. However, in a 
proper combustion model they have to be mutually re- 
lated since, when molecular processes are neglected, the 
flame front has to be solely fueled and carried by the tur- 
bulent dynamics of the reacting environment. To formu- 
late a correspondence between particle spread and flame 
progression, the expansion rate of the quadratic mean 




Figure 1: Illustration of the forward motion of the flame front 
from the initial position Co to Cf(t) and the expansion of 
the particle quadratic mean a after a forward plus backward 
random walk around the average position (x) = xo- 

of particle displacement r(r, t) = a(t + r)/a(t) is taken 
equal to the expansion rate of the flame front progression 
Af(t) = Cf{t) — Cq. In the mechanism here proposed, 
while the combustion evolves moving along the outward 
flame front normal, see Fig. [l] the expansion of a is 
statistically related to the random oscillation of a parti- 
cle moving forward and backward around its mean posi- 
tion. As a consequence, a half time step r is necessery to 
the forward-moving flame front to have the same expan- 
sion rate than an unburned particle oscillating around 
its mean position: r(r, t) = A/(t + r/2)/A/(t). Finally, 
noting that r(0,t) = 1 and performing the limit r — > 0, 
the joined process satisfies 
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dC 



f 



2 da 

C f (t) -C dt ~ alf)~dt 
Let us introduce the function 



Or 



T=0 



V{t) -2Ht 



I B L (r)dr, 
Jo 



(9) 



(10) 



where B^t) is the Lagrangian velocity autocorrelation 
function. It is worth remarking that definition ( 10 ) is in 

agreement with the exact Taylor formula cr 2 (t) — 2 J (t— 
t)Bl(t) dr that includes all turbulent dispersion regimes 
from the ballistic to the diffusive one passing through the 
inertial range. Using definitions |5]) and (10), identity ^ 
can be written in terms of U and T> and yields 



W) £v(0d£ 



(ii) 



It follows that f [V(t)U{£) - l4(t)V(£)]d£ = for any 
t > 0. From the monotonicity of both 14 and T>, since 
they are non-negative functions and 14(0) — T>(0) = 0, 
the following equality holds 



p(0 = 
w(0 u{t) 



< C < 



(12) 



so that, within the integration interval of £, the ratio 
T>(£)/14(£) must be constant. Moreover, for t — > oo both 



3 



V(t) and U{t) are bounded, V{oo) = V eq and U(oo) 
U eq ; therefore, such a constant is equal to 

V(t) V e 



U{t) u e 



A, t>0. 



(13) 



the 



Identity ([13j) states that, by using definition (10 1 
whole evolution of the combustion process is solely es- 
tablished by <r 2 (t), and it constitutes a new result in lit- 
erature [7]. Moreover, it determines not only the relation 
between the combustion drift and the background tur- 
bulent dispersion, but also the temporal evolution of the 
flame front location Cf{t) as it follows 



C f (t) = Co 



2A 



(14) 



A similar result was already sketched by Biagioli [TU] in 
the study of strongly swirled flows, but valid only for 
asymptotic long times. 

To validate the goodness of the physical argument that 
brought to expression (14 1, an experimental result dis- 
cussed in the literature is used where a and Cf were 
simultaneously measured. The data come from figure 
21 in Ref. [4 and they are attributed to unavailable 
measurements [11] . In the experimental setup, the gas 
mixture flows at the steady velocity V = 26.0 ms" 1 
and the data acquisition was performed at fixed dis- 
tances ie [mm] from an origin. The measurement loca- 
tions a; [mm] have been converted in elapsed times t[s] by 
t = x/V and Cf and a have been considered as mea- 
sured in a reference frame in translation with the flow. 
Figure [2] shows the fit between the measurements of a 
and Taylor formula, where the exponential autocorrela- 



tion function -Bi(i) 



)exp(— t/Tj,) is assumed with 



) T L = 11 X 10 a mm 2 s" 1 and T L = 2 x 10~ J s. 
Here T L = (u 12 )- 1 J Q B l {t) dr denotes the La grangian 
integral ti mesc ale. Substituting the resulting analytic 
a 2 in Eq. ( 14 ) , the flame front position Cf is correctly 
predicted. The experiment was performed twice with 
the same turbulence characteristics but with two differ- 
ent equivalence ratios F (that is the ratio of the fuel-to- 
oxidizer ratio to the stoichiometric fuel-to-oxidizer ratio): 
F = 0.68 and F = 0.56. 

Let us consider now the simple non-Markovian 
parabolic model: £ x [p] = T>{t)V 2 p. Equation Q be- 
comes the extension to all elapsed times of the familiar 
Zimont equation [3] , which was historically formulated in 
the asymptotic regime t ^> Tl with T>(t) tending to T> eq 
and U(t) to lA eq . A critical review about it can be found 
in [KB]. 

Despite the restriction n — 0, other extensions of Q 
to the initial regime t < Tl are currently applied in en- 
gineering applications to study transient and geometri- 
cal effects in the developing phase of the flame Q2] with 
practical fallouts in the design of spark-ignition engines 

With k — 0, the study of Zimont equation is reduced 
to a one-dimensional problem along the normal direction 
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Figure 2: Experimental validation of the propagation law 
(141. Lines are the plots of the analytic a and the predicted 
flame front positions Cf. The values A = 5.1mm and A = 
8.2 mm correspond to two different equivalence ratios F = 
0.68 and F = 0.56, respectively. 

to the flame front. The front speed Q becomes u(x, t) = 
— U(t) sign(dcfdx). Let the portion of space surrounded 
by the flame be fi(i) = [CL(t),Cn(t)], where Cl and Cr 
are the flame front positions defined in ^ on the left and 
on the right of the ignition point, respectively. Then the 
exact solution is 



c(x, t) 



jErfc 


' x - C R {t)~ 


-Erfc 


~x-c L (ty 











(15) 

By setting x — z — zq and Cn(t) — Zf(t) — zq, in the 
limit zq — > — oo the progress variable becomes c(z,t) = 
Erfc [(z — Zf(t))/(\/2 cr(i))] /2, that is in agreement with 
several experimental results [H [B]. In particular, for- 
mula (14 1 can be plugged into solution |l5| of the ex- 
tended Zimont equation. With this model, now fully 
determined for all elapsed times by a 2 , it is possible 
to perform a generally valid analysis on the flame en- 
hancement and quenching and bring out the different 
regimes of the combustion process. The two examples 
of Fig. [3] display the evolution of c(x, t) with the ini- 
tial fully burned zone bounded by different flame front 
conditions: £0 = A/5 and Co = 5 A. For the plots, 
we have set (u' 2 ) = 1, Tl = 1, A = 0.1. Observe in 
Fig. [3^1 that, after the initial ignition instant, turbu- 
lent diffusion mixes together products and reactants so 
that the progress variable takes a bell-shape profile and 
in the central zone c{x, t) is drained by the dispersion 
of the resulting product particles. If c(x, t) becomes less 
than a threshold value, quenching can be assumed. The 
larger is the ignition region the weaker is the draining, 
and, as illustrated in Fig. [3]d, no-draining occurs for an 
initial value Co larger than a critical lenghtscale Cq c . 
The lenghtscale £ c can be determined by the budget 
between the spread of the particle distribution, driven 
by T>(t), and the expansion of f2(i), driven by U(t). In- 
troducing a pseudo diffusion coefficient for the combus- 
tion U(t)Co, the budget emerges to be given by the ratio 
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Figure 3: Evolution of c(x, t) starting from a fully burned 
zone: a) Co = A/5 and t = 0, 0.01T L , 0.03Tl, 0.05T L , 0.07T l ; 
b) C = 5 A and t = O.OITl, O.ITl, 0.3T i; 0.5Tl, 0.7Tl. In 
a), the dashed line indicates the quenching threshold c = 0.5. 

U(t)Co/T>(t) = Cq/X. It emerges from this analysis that 
£oc = A = V eq /U eq = yJV eq T, where T = V eq /U% q is a 
determination of the characteristic reaction time. When 
£o « A, the flame rapidly quenches after the initial igni- 
tion, see Fig. [3^; when Co >• A, the flame is immediately 
able to sustain itself, see Fig. [3]d, and the combustion 
propagates. This dependence on the size of Co is the same 
as the one theoretically estabilished for the Kolmogorov- 
Petrovskii-Piskunov models [H]. In the intermediate 
case Co ~ A, when turbulence is not strong enough to 
drain c(x,t) below the quenching threshold, c(x, t), after 
an initial fall, is refilled and resustained by the combus- 
tion. This last behavior identifies the regime in which 
the process is dependent on the initial condition. Ac- 
tually, considering the right semiaxis on the flame front 
location x = Cn(t), the two contributions of the Erfc- 
functions in (15) become 1 and Erfc{v2 [Co/a+a/(2\)]}, 
respectively. Since a(t) is monotonic and increasing 
with <t(0) = 0, there exist two timescales Tq and t* 



defined by cr(r ) = Co and ct(t*) = 2A, respectively, 
which can be estimated from the scaling laws of c 2 (i). 
For short times, t <C T L , when a 2 ~ (u' 2 )t 2 , it turns 
out that To = £o/(w' 2 ) 1 / 2 . After this transient regime 
which is dependent on the initial condition, the process 
tends asymptotically to be self-similar. The timescale 
r* turns out to be r* - 2((?/ 2 ) /U^ 2 T L , if r, < T L , 
or = 2\/U eq , if r* 3> Tl. This means that when 
t < t q < then Cojo > cr/(2A) while when t > t* 
then £o/c < °7(2A). So when t < the second Erfc- 
function is strongly variable, while when t > r„ its argu- 
ment begins to grow and the Erfc-function tends rapidly 
to zero. Finally, when t < r* the average progress vari- 
able profile is not self-similar, because each Erfc-function 
in ( 15 1 has its own self-similarity variable, while for t > 



it is asymptotically self-similar with respect to the self- 
similarity variable (x — Cn(t))/a(t). Then the timescale 
r* separates the initial not self-similar transient regime 
from the asymptotically self-similar regime. For long 
elapsed times, in the reference frame moving with the 
flame front, the smoothing effect of the turbulence spread 
decreases and asymptotically vanishes, making the pro- 
file of c{x, t) steeper. The same argouments above can 
be applied to the left semiaxis. In addition, it is possi- 
ble to expand the present formulation including the ad- 
verting mean velocity by following the same arguments 
formulated in 15 where T> eq is changed in an effective 
diffusion coefficient with a nontrivial dependence on the 
mean velocity field. Then, by replacing ^ and ( 10 ) with 
a Lagrangian dispersion model for shear or cellular flow, 
the dependence of Cq c on the mean velocity may be es- 
timated as in [TBJ H3 ■ 
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